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Abstract 

Via exploitation of surface plasmon polaritons (SPPs), columnar thin films (CTFs) are attractive potential 
platforms for optical sensing as their relative permittivity dyadic and porosity can be tailored to order. 
Nanoscale model parameters of a CTF were determined from its measured relative permittivity dyadic, after 
inverting the Bruggeman homogenization formalism. These model parameters were then used to determine 
the relative permittivity dyadic of a fluid-infiltrated CTF. Two boundary-value problems were next solved: 
the first relating to SPP-wave propagation guided by the planar interface of a semi-infinitcly thick metal 
and a semi-infinitely thick CTF, and the second to the plane-wave response of the planar interface of a 
finitely thick metallic layer and a CTF in a modified Kretschmann configuration. Numerical studies revealed 
that SPP waves propagate at a lower phase speed and with a shorter propagation length, if the fluid has a 
larger refractive index. Furthermore, the angle of incidence required to excite an SPP wave in a modified 
Kretschmann configuration increases as the refractive index of the fluid increases. 
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1 Introduction 

Research on optical sensors exploiting surface plasmon polaritons (SPPs) has recorded an explosive growth 
[1, 2, 3] and commercial success [4, 5] during the last two decades. SPPs are virtual particles [6] that can 
propagate guided by the planar interface of a metal and a dielectric material, the classical analog of a train 
of SPPs being an SPP wave. 

Although several different sensing configurations have been investigated to excite SPP waves, perhaps the 
most commonplace is the Kretschmann configuration [7]. Both the metal and its dielectric partner in this 
configuration are layers of finite thickness. The metal film's thickness is ~ fOO nm, whereas the dielectric 
layer has to be much thicker. On the other side of the metal film is a dielectric coupling material (in the form 
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of a prism), which is optically denser than the dielectric material. Quasi- monochromatic light is launched 
at an angle to the thickness direction in the coupling material towards the metal film. The fraction of 
illuminating light that is neither reflected nor transmitted is absorbed. As the angle of incidence increases 
from towards 7r/2, a sharp peak in absorbance, accompanied by minuscule reflectance and transmittance, 
indicates the excitation of an SPP wave. This sharp peak occurs only for p-polarized light, and its angular 
location shifts when the refractive index of the dielectric partnering layer is altered by the incorporation of 
an analyte. 

The analyte can be incorporated in a variety of ways. One is to use a liquid such as water as the dielectric 
partner. Either the dispersal of the analyte in that liquid or the binding of the analyte to recognition 
molecules attached to the metal film can lead to an angular shift [2, 3]. Another way is to use a highly 
porous, thin metal film [8, 9], so that it can be infiltrated by an analyte-containing liquid [10, 11]. 

The dielectric partnering layer need not have an effectively isotropic constitution for the excitation of SPP 
waves. Uniaxial dielectric materials have been extensively considered for the role of the dielectric partnering 
material, both theoretically and experimentally [12, 13, 14, 15, 16]. Another class of candidates for this role 
is that of columnar thin films (CTFs), which are effectively biaxial dielectric materials [17]. 

A columnar thin film is an assembly of parallel nanowires, grown on a planar substrate by physical vapor 
deposition [17, 18]. As schematically illustrated in Fig. 1, the column inclination angle x € (0, 7r/2] is equal to 
or greater than the vapor incidence angle Xv G (0, 7r/2]. By control of the deposition-process parameters — 
the vapor incidence angle in particular — and the material(s) evaporated, both the optical properties and 
the porosity of CTFs may be tailored to order, thereby controlling the excitation conditions for SPP waves 
[19, 20]. This renders CTFs attractive as potential platforms for optical sensing of chemical and/or biological 
species that may infiltrate the void regions of the CTF [21, 22]. 

Our aim for the study reported here was to determine the sensitivity of SPP-wave excitation at the 
planar interface of a metal and a CTF to changes in the refractive index of the fiuid infiltrating the void 
regions of the CTF. Employing the relative permittivity dyadic of an uninfiltrated CTF, this study is largely 
numerical. It is primarily based on theoretical results established in two recent papers: one on the theory 
of exciting an SPP wave that propagates parallel to the morphologically significant plane of the CTF [19], 
the second on inverting the Bruggeman homogenization formalism with data on the uninfiltrated CTF [23] 
for use in a forward Bruggeman formalism [25] to compute the effective relative permittivity dyadic of the 
infiltrated CTF. 

The plan of this paper is as follows: Section 2 succinctly describes the inverse and forward Bruggeman 
homogenization formalisms for the uninfiltrated and the infiltrated CTF. Two boundary-value problems are 
formulated and solved next: 

(i) The first is a canonical problem relating to SPP-wavc propagation guided by the planar interface of a 
semi-infinitcly thick metal and a semi-infinitely thick CTF (Sec. 3), and 

(ii) the second is a realistic problem involving the plane-wave response of the planar interface of a finitely 
thick metallic layer and a CTF in a modified Krctschmann configuration (Sec. 4). 

Detailed descriptions of theoretical formulations are available in the predecessor papers [19, 23]. In the 
notation adopted, 3x3 dyadics are double underlined while 3-vectors are in boldface; unit vectors along the 
Cartesian axes are denoted u^,, and u^. The operators Re [•] and Im [•] deliver the real and imaginary parts 
of complex- valued quantities; and i = The free-space wavenumber is fco and the intrinsic impedance 

of free space is t^o- An exp(— iw^) time-dependence is implicit, with uj as the angular frequency and t as 
time. 
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2 Homogenization studies on CTFs 



2.1 Preliminaries 



Let us begin with an uninfiltrated CTF characterized at length scales much greater than the nanoscale by 
the relative permittivity dyadic 



efci u^u^ + eel UfcUfe. 



(1) 



Herein the normal, tangential and binomial basis vectors arc specified in terms of the column inclination 
angle x G (0, 7r/2] per 

u„ = -u^; sinx + cosx 



^ Uj; cos X + sm X 
Ub = -u„ 



(2) 



The CTF is grown on a planar substrate, which we take to lie parallel to z = 0, through the deposition of 
an evaporated bulk material. The vapor incidence angle Xv is less than or equal to the column inclination 
angle X: as schematically illustrated in Fig. 1. The deposited material is assumed to be an isotropic dielectric 
material with refractive index n^. The morphologically significant plane of the CTF is the xz plane [18, 
Chap. 7]. 

For the purposes of generating the numerical results presented in the following sections, we selected CTFs 
made from evaporating patinal® titanium oxide. The relative permittivity parameters of these CTFs are 



Cal 



ebi 



eel 



1.0443 + 2.7394 



1.6765 + 1.5649 



1.3586 + 2.1109 



2x. 

TT 

2x. 

TT 

2x. 

TT 



1.3697 



0.7825 



1.0554 



2x« 



2x. 



(3) 



ith 



tanx = 2.8818 tanx-u 



(4) 



as determined by experimental measurements at a free-space wavelength of 633 nm [24]. 

Each column of a CTF can be modeled as an assembly of elongated ellipsoidal inclusions, strung together 
end-to-end [18, Chap. 6]. The inclusions have identical orientations and shapes. The inclusion surface is 
characterized by the dyadic 

u„u„ +7^u^u^ + 7bUbUb. (5) 

The highly aciculate nature of the columnar morphology suggests that the shape parameters 7b ^ 1 and 
^ \. In fact, we fix 7t- = 15 since increasing 7t- beyond 10 does not bring about significant effects for 
slender inclusions [25]. The CTF's columns occupy only a proportion / G (0, 1) of the total CTF volume. 
The volume fraction of the CTF not occupied by columns is 1 — /. 



2.2 Inverse Bruggeman formalism — uninfiltrated CTF 

While the relative-permittivity parameters {e^i, ebi, eei} may be straightforwardly measured, the same is 
not true for the nanoscale model parameters {n.^,/, 7b}- The refractive index of the bulk material is known 
prior to evaporation, but the refractive index of the deposited material may be quite different, depending 
on the precise details of the deposition environment [26, 27, 28]. Also, the accurate determinations of 7b 
from scanning-electron micrographs [17], and / by gas-adsorption techniques [29, 30, 31] or mass density 
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measurements [26], are known to be problematic. Therefore, our first goal in this section is to establish the 
CTF nanoscale parameters {n^, /, 7^}, from a knowledge of the eigenvalues {^ai, £61, Cci} of e . The void 

— ct J ± 

regions of the uninfiltratcd CTF are taken to be vacuous. 

The computation of {ris,/, 7b} can be undertaken using the process of inverse homogcnization. The 
widely used Bruggcman homogcnization formalism [32, 33] provides an suitable framework for this process. 
Since we have recently described this inverse homogcnization process in detail elsewhere [23], here we need 
only present the results. 

For the titanium-oxide CTF described by the constitutive parameters (3) and (4), the corresponding 
values of {ug, f, jb} are presented in Table 1 for Xv S {15°, 30°, 60°, 90°}. The calculated data are in accord 
with previous observations: 

• As Xv increases towards tt/2, the biaxiality of a CTF reduces towards uniaxiality [17], which is reflected 
by the tendency of 7^ in Table 1 to reduce to almost unity. 

• CTFs deposited more obliquely tend to have lower mass density [18, Chap. 5], [34]. This is confirmed 
Table 1 by the monotonic growth in the volume fraction / as Xv increases towards tt/2. 

We also noticed a steady decrease in Ug towards the bulk refractive index of titanium oxide, as Xv it/2. 
This trend does depend on the material deposited, and can even be turned into a steady increase [23]. It 
suggests the complexity of the roles played by the nanoscale parameters {n^, /, 7f,} in the emergence of the 
relative-permittivity parameters {cai, e^i, Cci}- 

2.3 Forward homogenization — infiltrated CTF 

In an optical sensor exploiting the SPP-wave phenomenon, the void regions of the CTF have to be filled 
with a fluid of refractive index n^. As a result, e comprising the eigenvalues {fai, eti, Cci} changes to 

— ct J i 

e comprising new eigenvalues {£02, £62, ec2}; the column inclination angle x remains unchanged. From 

— ctj2 

knowledge of the nanoscale model parameters {?t-s,/, 7b} as well as of {ni,^r}, the Bruggeman homogcniza- 
tion formalism can be applied in its usual forward sense to determine {ea2, £62, £02}- A description of the 
forward Bruggeman formalism as applied to a CTF is available elsewhere [25] . 

After using the nanoscale model parameters presented in Table 1, the corresponding values of {ea2, £62, £c2} 
were plotted against £ (1.0, 1.5) in Fig. 2 for Xv <= {15°, 30°, 60°, 90°}. We see that each of €02, £62 and 
ec2 increases approximately linearly as ng increases. Also, the rate of increase is greater for smaller values of 
Xv , because the CTF is then more porous. 

3 Canonical boundary— value problem 

Let us now investigate the wavcnumbcrs of SPP waves guided by the planar interface of two half-spaces, one 
filled with a metal and the other occupied by a CTF. Since the underlying theory for this canonical problem 
is described elsewhere [19], here we present only the key theoretical results as a precursor to our numerical 
results. 

Suppose that the half-space z > is occupied by a fluid-infiltrated CTF while a metal of relative 
permittivity em occupies the half-space z < 0. The SPP wave under consideration is p-polarized and 
propagates in the xz plane. In the metal half-space, the electromagnetic phasors are 

E(r) = Am ( - — ) exp [ikn {ax - iq,nz)] 

\ qrn J 

H(r) = Am—^Viy exp [iko {ax - iqmz)] 

Vaqrn 

where A,„ is the complex-valued amplitude, qm = \/ a^ ^ e,„ and akgu^ represents the wave vector of the 
SPP wave. By choosing Re [qm] > 0, we ensure that this wave decays away from the interface z = 0. The 
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corresponding electromagnetic field phasors in the fluid-infiltrated-CTF half-space are 
iaqc - (ea2 - £62) sinx cosx 



E(r) = A, 



H(r) = A 



o-^ - (ea2 cos2 X + £62 sin^ x) 



o- (ea2 - £62) sin x cos x - iqc (ea2 cos^ x + £^2 sin^ x 



(£a2 COs2x + £62 sin^ x)] 



exp [i/cp (ca: + iqcz)] 

Uy exp [i/cg (fcc + igc-z) 



z > 0, 



where is the complex-valued amplitude and the quadratic equation [19] 

ql (£a2 COS^ X + £62 Sin^ x) + 2i0-9c (£a2 - £62) sin X cos X - (£a2 Sin^ X + £62 COS^ x) + €a2£62 = 

yields qc- The choice Re [qc] > ensures that this wave decays away from the interface. 
Upon applying the standard boundary conditions at z = 0, the dispersion relation [19] 



Emcr^ + iq,n<y (£a2 - £62) sinx cosx + [qniQc 



£m) (£a2 cos^ X + £62 sin^ x) 







(7) 
(8) 

(9) 



emerges, from which the relative wavenumber cr may be extracted by numerical means. 

For illustrative purposes, let us choose the metal occupying z < to be bulk aluminum, for which 
Cm = —56 + 2\i at a free-space wavelength of 633 nm. In Fig. 3, the real and imaginary parts of cr are 
plotted against ni for Xv G {15°, 30°, 60°, 90°}. The real part of cr increases approximately linearly as ni 
increases, and it also increases as Xv increases. From this we infer that the phase speed of the SPP wave 
decreases as increases and as Xv increases. Similarly, the imaginary part of a also increases as both ni 
and Xv increase, and from this we infer that the attenuation of the SPP wave increases as m increases and 
as Xv increases. 



4 Modified Kretschmann configuration 

We now turn to a realistic setup for launching SPP waves along the planar interface of a metal film and a 
CTF of finite thickness. We examine a modification [19] of the Kretschmann configuration [7] wherein the 
regions z < and z > is are occupied by homogeneous, isotropic, nondissipative, dielectric materials with 
relative permittivity scalars and ei = nj, respectively. The fluid-infiltrated CTF occupies the laminar 
region ^ z < Le, while the laminar region < z < Lm is occupied by a metal with relative permittivity 
£,„. A schematic representation of this modified Kretschmann configuration is provided in Fig. 4. 

Let us consider a p-polarized plane wave in the region z < 0, propagating towards the metal-coated CTF 
at an angle 6inc € (0,7r/2) to the z axis, as described by the electromagnetic field phasors 



( — Ua- cos t 



= — u„ exp [i {kx + JTd z cos 



exp [1(^2; + 

*inc)] 



erf z cos 



^inc) ] 



z < 0, 



with 

K = K\f^i sin6'i„c- 

The corresponding reflected and transmitted electromagnetic field phasors are given as 

Ere/ = {Vix COS Q^c + sin fl^nc) CXp \i [kx - y/e^ Z cos Oinc 



H 



re/ 



Uy exp \i {kx - sJTdZCO^Qinc)\ 



and 



Efr = ( — Uj- COS Bi, 



'?0 



U2 sin Otr) exp \i {kx + y/ee (z — Ls) cos ( 
Uj, exp [i {kx + y/ee (z — Ls) cosOfr)] 



z < 0, 



Z> Ly 



(10) 



(11) 



(12) 



(13) 
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respectively, where 




££ sin 9tr = (J. 



(14) 



The complex-valued reflection and transmission coefficients, namely Vp and ip, are determined by solving 
the related boundary-value problem. For full details, the reader is referred elsewhere [18, 19]. 
In studying surface waves at the mctal-CTF interface, a key parameter is the absorbance 



A characteristic signature of SPP-wave excitation is a peak in graph of Ap versus 0i„c: when the angle of 
incidence exceeds the critical angle in the absence of the metal film. 

For our numerical study, we chose L„i = 10 nm and Lx; = L^n + 1000 nm. The existence of a critical angle 
for total reflection in the absence of the metal film follows by setting the relative permittivity = 2.6^ = 6.76, 
in the neighborhood of that which can be delivered by a rutilc prism. In contrast, the relative permittivity 
ti was kept variable. As in Sec. 3, aluminum was chosen as the metal. 

In Fig. 5, graphs of Ap versus Oinc arc presented for ni G {1.0,1.25,1.5} with Xv — 30°. Also plotted 
versus Oinc is the quantity | p, computed when Lm = 0; these plots exhibit an abrupt step from to 1, 
which occurs at the critical value of Oinc for total reflection. The rightmost peak in the graphs Ap, which 
occurs beyond the critical angle, indicates the excitation of SPP waves. 

The fairly sharp spikes in the graphs of Ap, to the left of the rightmost peak, indicate bulk guided 
modes [35]. The Oinc values which correspond to these modes depend upon the thickness of the CTF. To 
demonstrate this, the plots of Fig. 5 are reproduced in Fig. 6 but with = i,„ + 1500 nm. From a 
comparison of Figs. 5 and 6, it is clear that bulk guided modes for = L,„ + 1500 nm arise at values of 
Oinc which are different to those for Lx; = + 1000 nm. Importantly, the values of Oinc which correspond 
to the excitation of SPPs are the same for — L„i + 1500 nm and is — Lm + 1000 nm. 

The value of Oinc which corresponds to the rightmost peak in the graph of Ap — let us call this value 
0\nc{''^i) a specified value of ni — is clearly sensitive to both ni and Xv- The value O^nd'^'i) increases as 
rii increases and as Xv increases. In order to quantify the sensitivity of Ol^^m), we introduce the figure of 
merit 



for ni E (1.0,1.5). Graphs of p versus ng are provided in Fig. 7 for Xv S {15°, 30°, 60°, 90°} and is = 
Lm + 1000 nm. These clearly indicate that the metal-coated CTF can function as an optical sensor for an 
analyte dispersed uniformly in a solution. Furthermore, we observe that 6'f„c("'^) niost sensitive to changes 
in ni when Xv is small and ng is large. 

5 Closing remarks 

Our numerical study has demonstrated that the excitation of an SPP wave guided by a planar interface of 
a metal films and a CTF is highly sensitive to the refractive index of a fluid infiltrating the CTF, as well 
as the morphology of the CTF itself. In particular, as the refractive index of the fluid increases, the phase 
speed of the SPP wave decreases and its degree of attenuation increases. Furthermore, the angle of incidence 
required to excite the SPP wave in a modifled Kretschmann configuration increases as the refractive index 
of the fluid increases. This sensitivity bodes well for the implementation of fluid-infiltrated CTFs, as well 
as fluid-infiltrated sculptured thin films [18] more generally, as SPP-bascd optical sensors. 
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P = 



ne - 1.0 



(16) 
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Xv 


lb 


./ 


Us 


15° 


2.2793 


0.3614 


3.2510 


30° 


1.8381 


0.5039 


3.0517 


60° 


1.4054 


0.6956 


2.9105 


90° 


1.0020 


0.7859 


2.8228 



Table 1: The dimcnsionlcss quantities 7f,, / and computed using the inverse Bruggcman homogenization 
formalism for a titanium-oxide CTF with Xv = 15°, 30°, 60° and 90°. 
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Figure 2: The relative permittivity-parameters 6^2, £62 and ec2 of the fluid-infiltrated CTF, as computed 
using the forward Bruggeman homgenization formalism, plotted against the refractive index of the fluid 
infiltrating the void regions of the CTF, for Xv — 15° (red, thick solid curve), 30° (green, dashed curve), 60° 
(blue, broken dashed curve) and 90° (black, thin solid curve). 
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Figure 3: The real and imaginary parts of a plotted against the refractive index of the fluid infiltrating 
the void regions of the CTF for Xv = 15° (red, thick solid curve), 30° (green, dashed curve), 60° (blue, 
broken dashed curve) and 90° (black, thin solid curve). 
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Figure 5: The absorbance Ap (red, solid curve) plotted against 6inc (in degree), when = 9, £,„ = — 56 + 21z, 
Lm = 10 nm, Ls = Lm + 1000 nm , Xv = 30° and ni £ {1.0, 1.25, 1.5}. Also plotted is the quantity |rpp 
(blue, dashed curve), calculated when i,„ = 0. 
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Figure 7: The figure of merit p (in degree) plotted against ni £ (1.0,1.5) for Ly, ~ L,n + 1000 nm with 
Xv = 15° (red, thick soHd curve), 30° (green, dashed curve), 60° (blue, broken dashed curve) and 90° (black, 
thin solid curve). 



16 



